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In the recent years, the interest in Stirling engine has remarkably increased due to its ability to use any 
heat source from outside including solar energy, fossil fuels and biomass. A large number of studies have 
been done on Stirling cycle analysis. In the present study, a mathematical model based on thermody¬ 
namic analysis of Stirling engine considering regenerative losses and internal irreversibilities has been 
developed. Power output, thermal efficiency and the cycle irreversibility parameter of Stirling engine 
are optimized simultaneously using Particle Swarm Optimization (PSO) algorithm, which is more 
effective than traditional genetic algorithms. In this optimization problem, some important parameters 
of Stirling engine are considered as decision variables, such as temperatures of the working fluid both 
in the high temperature isothermal process and in the low temperature isothermal process, dead volume 
ratios of each heat exchanger, volumes of each working spaces, effectiveness of the regenerator, and the 
system charge pressure. The Pareto optimal frontier is obtained and the final design solution has been 
selected by Linear Programming Technique for Multidimensional Analysis of Preference (LINMAP). 
Results show that the proposed multi-objective optimization approach can significantly outperform 
traditional single objective approaches. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Due to the dramatic consumption of fossil fuels, remarkable 
attention is drawn to the renewable energy and energy-efficient 
conversion systems. The world needs a clean energy revolution 
in order to break dependence on fossil fuels. Such a revolution 
would enhance globe energy security, promote enduring economic 
growth and tackle environmental challenges [1 ]. In this regard, the 
Stirling engine is one of the most promising sustainable energy 
technologies in recent years [2-4]. 

The Stirling engine is a simple type of external combustion 
engine which operates over a closed, regenerative thermodynamic 
cycle with the ability to use a wide variety of energy sources such 
as solar energy, geothermal and industrial waste heat. It plays an 
important role on environmental protection and also has signifi¬ 
cance for moderating the pressure on fossil fuels supplies in the 
world. With ideal regeneration and isothermal processes, the 
Stirling engine can theoretically convert heat into mechanical work 
at Carnot efficiency. The thermal limit for the operation of Stirling 
engine depends on the performance of heat resistant materials 
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used for its construction. In most instances, Stirling engines oper¬ 
ate with a heat source and heat sink temperature of 923 and 
338 K, respectively [5 . The engine thermal efficiency varies from 
about 30% to 40% resulting from a typical temperature in range 
of 923-1073 K, and normal operating speed range from 2000 to 
4000 rpm [6-10]. 

The invention of the Stirling engine preceded modern thermo¬ 
dynamic theory to a degree exemplified by the fact that it was 
not until fifty years after its invention that a reasonable analysis 
of the cycle was published. This analysis was made by Gustav 
Schmidt in 1871 [11] and has been reported in a more comprehen¬ 
sive form in the book of Israel Urieli 12]. The analysis, however, is 
limited in that it assumes isothermal working spaces and ideal 
heat exchangers, and thus it calculates the efficiency as being only 
temperature-dependent (i.e., Carnot efficiency). In 1960, Finkel- 
stein presented an ideal analysis, in which he assumed adiabatic 
working spaces but maintained the heat exchangers as ideal [13]. 
The effects of heat transfer and imperfect regeneration on the per¬ 
formance of the irreversible Stirling engine was investigated by 
Wu et al. [14 . As the efforts to reduce certain loss mechanisms 
can tend to increase the detrimental impact of others, complemen¬ 
tary approaches are needed. It is noted that a practical Stirling 
cycle always departs significantly from the ideal cycle, mainly 
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Nomenclature 



C v 

specific heat at constant volume (J/kg I<) 

t 3 

working fluid temperature in the hot space (I<) 

e 

regenerator effectiveness 

T 3 ' 

working fluid temperature at the regenerator 

Es 

Stirling engine thermal efficiency 


outlet (I<) 

I< 

factor defined by Eq. (9) 

n 

working fluid temperature in the cold space (K) 

k 

specific heat ratio 

Tv 

working fluid temperature at the regenerator inlet (I<) 

1<DH 

hot space dead volume ratio 

t r 

effective working fluid temperature in regenerator dead 

l<DR 

regenerator dead volume ratio 


space (I<) 

kDC 

cold space dead volume ratio 

Vdh 

hot space dead volume (m 3 ) 

m 

total working fluid mass contained in the engine (kg) 

Vdr 

regenerator dead volume (m 3 ) 

P 

charge pressure (N/m 2 ) 

^DC 

cold space dead volume (m 3 ) 

Qin 

total heat added from an external heat source (J) 

V D 

total dead volume (m 3 ) 

Qout 

total heat rejected to an external heat sink (J) 

v E 

expansion volume (m 3 ) 

R 

gas constant (J/kg I<) 

v c 

compression volume (m 3 ) 

Rs 

irreversibility parameter 

W ne t 

engine net work (J) 


due to the performance of the heat exchanger and the regenerator. 
A detailed analysis with practical values and examples of the 
engine were given by Organ [15]. The similarity and scaling 
approach developed by Organ provides a simpler way of creating 
the preliminary design for the prototype engine. In 2000, Kaushik 
and Kumar presented an investigation of a finite-time thermody¬ 
namic analysis of a Stirling engine in which the power output 
and the corresponding thermal efficiency were maximized [16]. 
Petrescu and Costea [17 conducted a technique for calculating 
the efficiency and power of Stirling engines based on the First 
law of thermodynamics for processes with finite speed [17,18 . In 
2008, Timoumi and Tlili developed a second-order Stirling model 
which includes thermal losses and is applied for the optimization 
of GPU-3 Stirling engines [19]. Leonardo Scollo [20] presented an 
alpha-type Stirling engine by using a similar design method. A 
numerical model for beta-type Stirling engine with rhombic drive 
mechanism was developed by Cheng et al. [21]. The energy equa¬ 
tions of the control volumes in the working chambers and heat 
exchangers were derived and solved by taking into account the 
non-isothermal effects, the thermal resistance of the heater head 
and the effectiveness of the regenerative channel [21 . Formosa 
et al. [22] developed an analytical model taken into account the 
heat losses and irreversibilities on Stirling engines which could 
be used in a dynamical analysis for preliminary design [22]. In 
2011, a semi-analytical dynamic model of free piston Stirling 
engine was developed by Formosa [23]. The thermodynamic model 
is used to define the thermal variables which are used in dynamic 
model which evaluates the kinematic results. Recently, Aksoy and 
Cinar 24] conducted a theoretical investigation on dynamic and 
thermodynamic analysis of a beta-type Stirling engine with rhom¬ 
bic drive mechanism by implementing nodal analysis method [24]. 
Solmaz and Karabulut [25] presented a novel configuration of beta- 
type Stirling engine which was driven by a lever mechanism. The 
performance of the novel engine was compared with rhombic drive 
engine via nodal analysis [25]. Abdollahpour et al. [26] developed a 
thermodynamic model and an optical model of a parabolish dish 
collector for Stirling engine. The collector has the ability to be cou¬ 
pled with different type of Stirling engine and used in various fields 
like air conditioning and solar water heater in the following of gen¬ 
erate electricity from Stirling engine [26]. 

Although a large number of studies have been done on the opti¬ 
mization of Stirling cycle and Stirling engine, most of them are 
merely focused on single impact factor. Previous work has been 
done on dynamic and thermodynamic analysis of Stirling engine 
by the authors [27,28]. In this paper, we have studied three 
objective functions simultaneously: engine power output, thermal 
efficiency of Stirling engine, and the cycle irreversibility parameter 


which has previously been proved to be significant in predicting 
the performance of Stirling engine in the thermodynamic cycle 
[29]. Meanwhile, ten decision variables including temperatures of 
heat source and heat sink, volumes of expansion space and com¬ 
pression space, total dead volume, volume ratios of three heat 
exchangers, effectiveness of regenerator, and mean effective pres¬ 
sure are considered during the multi-objective optimization 
process. 

Problems with multiple objectives are present in a great variety 
of engineering optimization problems. It is quite common that in 
these problems there are several conflicting objectives to be opti¬ 
mized and it is difficult to identify what the best solution is. 
Despite the considerable diversity of techniques in Operation 
Research developed to tackle these problems, their intrinsic com¬ 
plexity calls for alternative approaches. Recently, Ahmadi et al. 
[30-35 conducted a series of works on multi-objective optimiza¬ 
tion method for designing a power Stirling engine by implement¬ 
ing the second-version Non-dominated Sorting Genetic Algorithm 
(NSGA-II) which is known as a tributary of evolutionary algorithm. 
In this study, a superior technique that has been adopted for deal¬ 
ing with multi-objective optimization problems is Particle Swarm 
Optimization (PSO) [36 , which is a relatively heuristic inspired 
by the choreography of a bird flock. This method was applied in 
the area of designing Stirling engine for the first time. The study 
developed by Coello [37] demonstrated that PSO is highly compet¬ 
itive and can be considered a viable alternative to solve multi¬ 
objective optimization problems. In addition, the exceptionally 
low computational times required by PSO approach make it a very 
promising approach to problems (e.g., engineering optimization) in 
which the computational cost is a vital issue. 

2. Stirling system description 

The ideal Stirling cycle satisfies the Carnot requirement of 
reversibility and can be described with reference to Fig. 1. One of 
the mechanical configurations for realizing the Stirling cycle is 
shown, being known as the Alpha arrangement. The regenerator 
generally comprises a matrix of fine wires, porous metal or some¬ 
times simply the metal wall surface enclosing an annular gap. 

The position of the pistons are shown at the four extreme state 
points of the cycle as seen in the pressure-volume and tempera¬ 
ture-entropy diagrams. Process 1-2 is the isothermal compression 
process during which the heat is removed from the engine at the 
cold sink temperature. Similarly, process 3-4 is the isothermal 
expansion process during which heat is added to the engine at 
the hot source temperature. Process 2-3 and 4-1 are the constant 
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Fig. 1. Illustration of the Stirling cycle system. 


volume displacement processes in which the working gas is passed 
through the regenerator. During the process 4-1 the working gas 
gives its heat up to the regenerator matrix, and the heat will be 
recovered from the regenerator matrix to the gas subsequently 
during the process 2-3. Thus, even though quite considerable 
quantities of heat are transferred during the displacement process, 
they are seen to be externally adiabatic, fulfilling the Carnot 
requirements for maximum attainable efficiency. 


In a similar way, the working fluid temperature at the regener¬ 
ator inlet can be expressed as: 

TV = T 3 + e(J\ - T 3 ) (4) 

For a real Stirling cycle, the accurate description of working 
fluid temperature for the regenerator is important. The effective 
temperature of the working fluid within the regenerator dead vol¬ 
ume can be determined using the arithmetic mean [38]: 


3. Analysis of Stirling cycle with internal irreversibility 

3.1. Stirling engine dead volumes 

Total dead volume includes the heater, the regenerator and the 
cooler dead volumes and is expressed as: 

= Vdh + V dr + Vdc = (kDH + k D R + I<DcWd (1) 

where k DHt k DR and k DC are the heater dead volume ratio, the regen¬ 
erator dead volume ratio and the cooler dead volume ratio, 
respectively. 

3.2. Imperfect regenerator 


TV + TV 

Tr= ^~ 

Substituting Eqs. (3) and (4) into Eq. (5) gives: 

T h+T, 

Tr= ^2— 


( 5 ) 

( 6 ) 


3.3. State equation 

Assume that the expansion and compression volumes are, V E 
and V c respectively. Using the ideal gas law given by: 

pV = mRT (7) 


For an imperfect regenerator, the working fluid temperature at 
the outlet and inlet of the regenerator are T 3 - and TV, respectively. 
The regenerator effectiveness, e, is defined as follows [38]: 



The regenerator effectiveness e thus varies from 0 for no regen¬ 
eration to 1 for an ideal regenerator. The temperature of the work¬ 
ing fluid at the regenerator outlet can be expressed in terms of the 
regenerator effectiveness as: 

T 3 '=Ti+e(T 3 -Ti) (3) 


the equation of state for the isothermal compression process 1-2 
with dead volumes V DH , V DR and V DC can be expressed as follows: 


P 


mR 

Y_e_ i Vdh i Vdr i Vdc 
t 3 "t" r 3 "r" t r _r t i 



mR 


(g + * + 




where 


K 


Vdh V D r , 
+ + 


Vdc 

tt 



It should be noted that, for a given hot-side and cold-side working 
fluid temperature, the factor K is a function of the dead volumes. 
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3.4. Thermodynamic processes of Stirling cycle 

For the isothermal compression process 1-2, the heat rejected is 
given by: 


W ne t — Q. — Qin Qout 


pi 




T 


T 3 In 


i 


(w+w+/<r 3 y 

(V e + I<T 3 ) 


Ti In 


(W + W+KTi) 


(V £ + / 0 ’ 1 ) 


( 20 ) 


>v f 


0 . 1—2 = VVi_2 


'V E +Vc 


pdV c = mRT i In 


(Vjr + KTO 
(\/ £ + Vc + /<Ti)J 


( 10 ) 


It can be seen that the compression work depends only on the 
factor I< which is a function of the dead volumes. 

For the isochoric heating process 2-3, the heat added is given 


Q 2—3 = mCv(T 3 — T 2 ) = mCy(T 3 — T 1 ) 



The engine thermal efficiency can be determined from: 


w„ et _ ■ 

>3 In 

\V E +V c +KT 3 j 
( V e +I<T 3 ) 

- 7, In 

-(V E+ V c +l<T,j 
(V e +KT 1) 

I 

Qin 

T 3 In 

\V E +V c +KT 3 j 

( v e +i<j 3 ) 

(l-e)(r 3 -T 1 )'| 
“r k-\ J 



( 21 ) 


In order to describe quantitatively the effect of the internal dis¬ 
sipations of the working fluid on the performance of Stirling 
engine, the cycle irreversibility parameter, R s , is defined as: 


where C v is the specific heat at constant volume and is assumed to 
be constant. Without regeneration, this amount of heat is added 
from an external heat source, while for ideal regeneration this 
amount of heat is released from an ideal regenerator. 

The regeneration heat released from an imperfect regenerator 
during the process 2-3’ is: 

0 .2— 3 ' = tnCy(T 3 ' — T 2 ) = emCv(T 3 — T\) (12) 

The heat added from the external heat source during process 3’- 
3 is: 



(V D + Vp + i<t 3 \ / fv D + Vp + icr\\ 
V Vd + I<T 3 )/ V Vo + KT, ) 


S4 — S3 

Si — s 2 


( 22 ) 


It shows that the internal irreversibilities in a Stirling engine 
can be characterized by a single parameter representing the ratio 
of two entropy differences. This parameter has not been widely 
adopted in designing Stiring engine. In this study, R s is considered 
as one of the objectives to be maximized. It has been proved in [29] 
that Stirling engine with higher the cycle irreversibility can deliver 
more power and has higher efficiency. 


0 3 '- 3 = mC v (T 3 - TV) = (1 - e)mC v (J 3 - T a ) (13) 


It is clear that the heat input to the isochoric process depends only 
on the regenerator effectiveness. 

In the isothermal expansion process 3-4, the heat added to the 
cycle is given by: 


0 . 3 -4 



mRT 3 In 


\(V E + Vc + KT 3 )] 

(v E + tcr 3 ) 



It is obviously that expansion work depends on the dead volumes. 
For the isochoric cooling process 4-1, the heat rejected is: 

O 4—1 = TnCv(T\ — T 4 ) = —mCy(T 3 — T\) (15) 

Without regeneration, this amount of heat is rejected to an 
external heat sink, while for ideal regeneration this amount of heat 
is rejected to the ideal regenerator. 

For an imperfect regeneration, the heat rejected to the regener¬ 
ator is given by: 

Q. 4 _i' = mC v (T v - T 4 ) = -emC v (T 3 - Ti) (16) 

The heat rejected to an external heat sink during process l’-l is 
given by: 

Qi'_i = mC v (h - TV) = -(1 - e)mC v (J 3 - T t ) (17) 

It is evidenced that the heat transfer during the cooling process 
depends only on the regenerator effectiveness. 

Consequently, for an imperfect regeneration, the total heat 
added and heat rejected in Stirling cycle can be expressed as 
follows: 



Q.3— 3 ' + Q 3-4 

mC v |(l -e)(T 3 -T 1 ) + (fc-l)T 3 ln 

(18) 


( V e + V c + KT 3 ) 
(V E + KT 3 ) 


lout — Q. 1 - 1 ' + O 1-2 

= -mCvi (1 - e)(T 3 - h) + (k - 1)1, In 


(V e + Vc + I<T 1 ) 
(Ve + KT,) 


(19) 


The engine net work can be determined from: 


4. Multi-objective optimization via particle swarm optimization 
algorithm 

4.1. Multi-objective particle swarm optimization 

The optimization of Stirling engine is a typical multi-objective 
optimization problem. Lets c 91 n be an n-dimensional search space, 
and fix), i = 1 be k objective functions which are defined over 

5. Also let / be a vector function defined as, f(x) = [fj (x),/ 2 (x),..., 
fk(x)], and gfx) ^ 0, 1 = 1,..., m be m inequality constraints. Then 
the goal of multi-objective optimization is to find a solution, 
X* = (X\,X*2, ■ • • >X*)> that minimizes/(x). 

To solve this problem, an advanced technology in optimization 
is adopted, which is Multi-objective Particle Swarm Optimization 
(MOPSO) [37]. MOPSO is a general optimization algorithm, but it 
has never been applied in the research area of Stirling engine. In 
this study, MOPSO is proved to work extremely well in task of 
designing Stirling engine. 

In order to introduce MOPSO, a review of Particle Swarm Optimi¬ 
zation (PSO) algorithm is given. As a classic optimization technique, 
PSO was first proposed by Kennedy and Eberhart in 1995 [36 . It is a 
computational method that optimizes a problem by iteratively try¬ 
ing to improve a candidate solution with regard to a given measure 
of quality, e.g., an objective function in our case. This method is sim¬ 
ilar in some aspects to other evolutionary algorithms with the 
exception of that particles move rather than evolve through the 
search space. If enough particles are given, PSO guarantees to get 
the global optimum solution of the objective function. In practice, 
however, fewer particles are used to find a balance between accu¬ 
racy and efficiency, which usually leads to suboptimum solution. 

As compared to Genetic algorithms (GA), there are two main 
differences between PSO and GA (39 . Firstly, genetic algorithms 
rely on three mechanisms in their processing: parent representa¬ 
tion, selection of individuals and the fine tuning of their parame¬ 
ters. On the contrary, PSO only relies on two mechanisms 
because it does not adopt an explicit selection operator. However, 
the absence of a selection mechanism in PSO is compensated by 
the use of leaders to guide the search. A second distinction has to 
do with the way in which the individuals are manipulated. PSO 
uses an operator that sets the velocity of a particle to a particular 
direction. This can be seen as a directional mutation operator in 
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which the direction is defined by both the particle’s personal best 
and the global best of the swarm. If the direction of the personal 
best is similar to the direction of the global best, the angle of 
potential directions will be small, whereas a larger angle will pro¬ 
vide a larger range of exploration. In contrast, GA uses a mutation 
operator that can set an individual in any direction, although the 
relative probabilities for each direction may be different. 

In multi-objective optimization, the objectives to be optimized 
are normally in conflict to each other, which means that there is no 
single solution for these problems. Instead, MOPSO aims to find 
good tradeoff solutions that represent the best possible compro¬ 
mises among the objectives based on the concept of Pareto opti¬ 
mality. A point X* c Q is Pareto optimality if for every X c Q 
and / = {1, 2,..., k} either 

ViciO m =£(>?)) (23) 

or, there is at least one i c I such that 

m>M?) ( 24 ) 

It is defined that X* is Pareto optimal if there exists no feasible 
vector X which would decrease some criterion without causing a 
simultaneous increase in at least one other criterion. In MOPSO, 
all objective functions are evaluated for each particle. Pareto opti¬ 
mality can guide the particles to produce nondominated best posi¬ 
tions (often called leaders). Since there can be many nondominated 
solutions in the neighborhood of particle, the determination of 
leaders is not straightforward, but only one is usually selected to 
participate in the velocity update in PSO. The whole scheme of 
MOPSO can be outline in Fig. 2, in which, the external archive is 
used for storing the nondominated solutions discovered during 
search [40-42 . An external archive has bounded size, thereby 
making unavoidable the imposition of the rules regarding the 
replacement of existing solutions with new ones. 

4.2. Objective functions , decision variables and constraints 

In the present study, three key objective functions for optimiza¬ 
tion are the output power (W net ), the Stirling engine thermal effi¬ 
ciency (E s ) and the cycle irreversibility parameter (R s ), denoted 
by Eqs. (20)-(22), respectively. There are ten decision variables 
considered as follows: 



Initialize swarm, velocities, best 
positions, and external archive 
_ (empty) _ 

i 

Stopping criterion 
satisfied? 


I 


no 


For each particle, select a 
member of external archive 


yes 


I 

For each particle, update velocity 

and position 


I 

For each particle, evaluate new 
position position 


i. 

For each particle, update the 
best position & external archive 



Tj\ Working fluid temperature in the cold space 

T 3 : Working fluid temperature in the hot space 

V D : Total dead volume 

V E \ Expansion volume 

V c : Compression volume 

e: Effectiveness of regenerator 

k DH : Elot space dead volume ratio 

k DR : Regenerator dead volume ratio 

k DC : Cold space dead volume ratio 

p: System charge pressure. 


The objective functions have been solved by implementing the 
MOPSO algorithm described above with respect to the following 
constrains: 

288 < r, < 320 I< 

(25) 

800 < r 3 < 1300 K 

(26) 

50 < V E < 150 cm 2 

(27) 

50<V C < 150 cm 2 

(28) 

75 < V D < 180 cm 2 

(29) 


Fig. 2. The flow chart of MOPSO. 


0.4 < e < 0.9 

(30) 

0.1 ^ ^ 0.2 

(31) 

0.2 ^ k[) R ^ 0.4 

(32) 

0.4 < k DC ^ 0.7 

(33) 

koH + 1<DR + k D c = 1 

(34) 

1 ^ p ^ 5 Mpa 

(35) 


In this study, a MOPSO based Matlab program is developed to 
find the optimum design parameters of the Stirling engine. The 
population size of MOPSO algorithm options is set to be 250. 


5. Decision-making methods in multi-objective optimization 

In multi-objective optimization, a decision-making process for 
selection of the best optimal solution from all the feasible alterna¬ 
tives is required. There exist a number of several methods for deci- 
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Table 1 

Specifications of the Stirling engine. 


Work gas 

Nom. molar mass 

Gas constant 

Specific heat 


Specific heat ratio 


M (kg/kmol) 

R (kj/kg I<) 

c p 

C v 

ii 

G 1 

He 

4 

2.08 

5.19 

3.11 

1.67 


pareto 



0.55 1000 


Irreversibility parameter Rs 


Power outputWnet 


200 


Fig. 3. Pareto optimal frontier in objectives’ space. 


Pareto optimal frontier in objective space 



sion-making process in Decision problems. And these methods can 
be implemented in decision-making for selection of a final optimal 
solution from the Pareto frontier which is obtained by MOPSO. 
Linear Programming Technique for Multidimensional Analysis of 


Preference (LINMAP), which is one of the most recognized and 
popular type of decision-making processes, is based on pair-wise 
solutions’ comparisons obtained from decision-making processes 
and yields the best compromise solution which has the shortest 
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distance to the positive ideal solution [43,44 . Recently, Ahmadi 
et al. [31,32 observed that the LINMAP decision-making method 
could generate a final optimal solution with the least deviation 
index from the positive ideal solution. Hence, in the present study, 
LINMAP is implemented and the final optimal solution has been 
decided based on this method. 

6. Results and discussion 

The engine power output, thermal efficiency of Stirling engine 
and irreversibility parameter are maximized simultaneously by 


implementing multi-objective optimization method based on par¬ 
ticle swarm optimization algorithm. 

In this study, the optimization process is conducted with three 
objective functions which are expressed by Eqs. (20)-(22) and ten 
design parameters that are restrained by Eqs. (25)-(35). The deci¬ 
sion variables (design parameters) of optimization are temperature 
of working fluid in the hot space, temperature of working fluid in 
the cold space, total dead volume, expansion volume, compression 
volume, dead volume ratio of the hot space, dead volume ratio of 
the cold space, dead volume ratio of the regenerator, effectiveness 
of the regenerator, and the system charge pressure. In this case, 


Pareto optimal frontier in objective space 



Pareto optimal frontier in objective space 



Fig. 6. Pareto optimal frontier in objectives’ space ( E s - Rs ). 
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Table 2 

Final optimal solution obtained by LINMAP decision-making method. 
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Methods 

Decision variables 









Objective functions 



h (K) 

T 3 (K) 

V D (cm 3 ) 

14 (cm 3 ) 

V c (cm 3 ) 

e 

l<DH 

l<DR 

kDC 

P( Mpa) 

Power (W) 

E s 

Rs 

LINMAP 

297.9 

914.1 

143.3 

98.1 

75.9 

0.89 

0.38 

0.41 

0.21 

4.61 

732.3 

0.18 

0.67 

Results [28] 

314 

923 

136.2 

96.8 

93.1 

0.88 

0.4 

0.5 

0.1 

4.1 

633.6 

0.15 

- 


helium is chosen as the working fluid. In order to have consistency 
with previous works [28], specifications of the Stirling engine are 
summarized in able 1. 

The Pareto optimal frontier in the proposed objectives’ space 
obtained by multi-objective particle swarm optimization is illus¬ 
trated in Fig. 3. The final solution selected by LINMAP decision¬ 
maker is also indicated in Fig. 3. Moreover, a surface is fitted to 
the points obtained by multi-objective optimization as following: 

E s = 6.25 x 10 4 - 1.09 x 10 5 K s -2.69 x 10 5 W, let + 2.43 
x 10 5 J?s + 2.48 x 10 5 R s W„ et +4.03 x 10 5 W^ t - 2.77 
x 10 5 R| - 2.39 x 10 5 RsW, Kt - 1.38 x 10 5 R s W 2 net 
-2.05 x 10 5 W^ t (36) 

which can be applied to gain the irreversibility parameter. 

In Fig. 4 Pareto optimal frontier for objective functions of power 
output and thermal efficiency of the Stirling engine is depicted. As 
can be seen in Fig. 4, the power output of the Stirling engine 
decreases when the thermal efficiency increases. In addition, the 
fitted curve equation of Pareto optimal frontier is obtained as 
following: 

W net = 1.77 x 10 7 £j - 9.48 x 10 6 Es + 1.69 x 10 6 E S 

-98479.11 (37) 

Pareto optimal frontier for power output and irreversibility 
parameter is presented in Fig. 5. It can be seen that engine power 
output decreases with increasing the irreversibility parameter of 
the Stirling cycle. To assist the optimal design of the Stirling 
engine, the fitted curve, which is valid in the range of 
0.64 <R S < 0.75, is derived as following: 

W„ e[ = -6068.98Rj + 4886.83E S + 440.56 (38) 

Fig. 6 illustrates the Pareto optimal frontier for thermal effi¬ 
ciency and irreversibility parameter of the Stirling engine. It is 
clear that thermal efficiency increases to 0.198 when irreversibility 
parameter decreases to 0.742. The fitted curve equation of Pareto 
optimal frontier is expressed as following: 

£ s = -12.03E 7 + 17.31E S - 6.02 (39) 

The final optimal solution for design parameters and objective 
functions using LINMAP decision-making method is presented in 
Table 2 compared with the experimental data [28] in which the 
cycle irreversibility parameter R s is not considered. The improved 
design outperforms the previous work 28] by more than 15% 
and 20% on engine power output and thermal efficiency, respec¬ 
tively. It is obviously that the design proposal obtained from the 
proposed MOPSO algorithm leads to better performance of Stirling 
engine than that of the traditional design. 

7. Conclusions 

According to the previous study on the Stirling cycle and Stir¬ 
ling regenerator, the thermodynamic analysis of the Stirling engine 
proposed in this paper shows more accurate results compared with 
the others due to the consideration of the cycle irreversibility 
parameter. The engine power output, thermal efficiency and the 


cycle irreversibility parameter are optimized simultaneously for 
the multi-objective optimization while temperature of working 
fluid in the hot space (T-i), temperature of working fluid in the cold 
space (T 3 ), total dead volume (V D ), expansion volume (V E ), com¬ 
pression volume (Vc), dead volume ratio of the hot space (k DH ), 
dead volume ratio of the cold space (k DC ), dead volume ratio of 
the regenerator (k DR ), effectiveness of the regenerator (e), and the 
system charge pressure (p) are considered as design parameters. 
Multi-objective optimization is conducted by implementing parti¬ 
cle swarm optimization algorithm and the Pareto optimal frontier 
in objectives’ space is obtained. Moreover, the final solution is 
selected from Pareto optimal by LINMAP decision-maker. The indi¬ 
cated power output is 732.3 W. The thermal efficiency and the irre¬ 
versibility parameter of the Stirling engine is 0.18 and 0.67, 
respectively. The simulation results show the correctness of our 
physical model and mathematical model, which suggests that the 
proposed multi-objective multi-variable optimization method is 
capable in designing Stirling engine. The performance of Stirling 
engine designed by MOPSO method is 15% and 20% better than 
previous work based on traditional single-objective approach on 
power output and thermal efficiency, respectively. 
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